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ABSTRACT 

We present a first application of the subhalo abundance matching (SHAM) method to describe 
the redshift-space clustering of galaxies including the non-linear redshift-space distortion, i.e., 
the Fingers-of-God. We find that the standard SHAM connecting the luminosity of galaxies 
to the maximum circular velocity of subhalos well reproduces the luminosity dependence of 
redshift-space clustering of galaxies in the Sloan Digital Sky Survey in a wide range of scales 
from 0.3 to 40 ft,“^Mpc. The result indicates that the SHAM approach is very promising 
for establishing a theoretical model of redshift-space galaxy clustering without additional 
parameters. We also test color abundance matching using two different proxies for colors: 
subhalo age and local dark matter density following the method by Masaki et al. (2013b). 
Observed clustering of red galaxies exhibits much stronger Fingers-of-God effect than blue 
galaxies. We find that the subhalo age model describes the observed color-dependent redshift- 
space clustering much better than the local dark matter density model. The result infers that 
the age of subhalos is a key ingredient to determine the color of galaxies. 

Key words: galaxies: formation - galaxies: haloes - galaxies: statistics cosmology: observa¬ 
tions - cosmology: large-scale structure of Universe 


1 INTRODUCTION 

Establishing connection between the properties of galaxies and the 
underlying dark matter is crucial for both studies of galaxy evolu¬ 
tion and cosmology. Star formation histories of galaxies has been 
studied by associating galaxies with their host dark matter halos 
and their connection provides fundamental constraints on galaxy 
formation models (e.g., Conroy &. Wechsler 2009; Leauthaud et al. 
2010; Behroozi et al. 2013). Future galaxy surveys such as Prime 
Focus Spectrograph (PFS) (Takada et al. 2012), the Dark Energy 
Spectroscopic Instrument (DESI) (Levi et al. 2013), Euclid (Lau- 
reijs et al. 2011) and the Wide Field Infrared Survey Telescope 
(WFIRST) (Spergel et al. 2013) use both luminous red galaxies and 
emission line galaxies to trace the large-scale structure at z 2. A 
major uncertainty for the precision cosmology using galaxy surveys 
comes from the challenge of relating galaxies and dark matter. 

Subhalo abundance matching (SHAM) is a promising ap¬ 
proach to relate the properties of galaxies to dark matter subhalos 
(e.g., Kravtsov et al. 2004; Nagai & Kravtsov 2005; Conroy et al. 
2006). The simple abundance matching model by assigning lumi¬ 
nosity in the order of the maximum circular velocity of dark matter 
subhalos successfully reproduces the galaxy clustering at different 
redshifts (Conroy et al. 2006). Masaki et al. (2013a) also finds a 
good abundance matching between the progenitor halos of lumi¬ 
nous red galaxies (LRGs) and the massive halos at 2 ; ~ 2 and then 
explains the clustering properties of LRGs very well. There also has 
been recent attempts to relate the galaxy color to the subhalo prop¬ 


erties, i.e., color abundance matching. Galaxy color reflects the ac¬ 
tivity of on-going star formation: red galaxies consists of aged stars 
and the star formation is quenched, whereas blue galaxies are rela¬ 
tively young and their star formation is active. It is also known that 
redder galaxies live in denser environments via the measurement of 
galaxy clustering (Norberg et al. 2002; Zehavi et al. 2005; Coil et al. 
2008; Guo et al. 2014; Skibba et al. 2014) and also from the color- 
density relation (Balogh et al. 2004; Cooper et al. 2006; Blanton 
& Berlind 2007; Bamford et al. 2009). Masaki et al. (2013b) ex¬ 
tends a SHAM technique to explain color-dependent properties of 
galaxy clustering as well as galaxy-galaxy lensing using two prox¬ 
ies of color: one is the local dark matter density motivated by the 
environmental dependence of galaxy color; the other is the subhalo 
age reflecting the different aged population between red and blue 
galaxies. Hearin & Watson (2013) also perform color abundance 
matching by assigning the redshift Zstarve characterizing the epoch 
of star formation quenching to subhalos. 

So far, the projected correlation function along the line-of- 
sight has been commonly used for testing SHAM to avoid the ef¬ 
fect of redshift-space distortion (RSD) due to peculiar motion of 
galaxies. The velocity of galaxies within and outskirts of clusters 
is complicated and affected by different physics including the dy¬ 
namical friction, tidal stripping/disruption, merging, and ram pres¬ 
sure. The internal motion of galaxies elongate the RSD of galaxies 
along the line-of-sight direction, known as Fingers-of-God (FoG; 
Jackson 1972). The FoG effect is clearly different by colors: red 
galaxies show much stronger FoG effect than blue galaxies (Zehavi 
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et al. 2005; Coil et al. 2008). In the framework of the halo model, 
the line-of-sight elongation in the galaxy distribution emerges from 
the one-halo term, which is the contribution of central-satellite 
and satellite-satellite galaxy pairs in the same halos (Hikage & 
Yamamoto 2013). The observed color difference of FoG feature 
mainly reflects the difference of the satellite fraction, radial pro¬ 
files and their kinematics between red and blue galaxies. 

In this paper we extend the SHAM approach to redshift space 
and test if the model describes the luminosity and color depen¬ 
dence of redshift-space galaxy clustering for the first time. We 
characterize the anisotropy of redshift-space galaxy clustering by 
a multipole expansion. High-i multipole components such as hex- 
adecapole {I = 4) mainly generated from the FoG effect provide 
a useful probe of constraining the fraction and the internal veloc¬ 
ity dispersions of satellite galaxies (Hikage 2014). In our analy¬ 
sis, we focus on the non-linear scales from subMpc to ~ IMpc 
where the difference of the FoG effect can be clearly seen. Extend¬ 
ing SHAM to redshift space is an important step for establishing a 
theoretical model of the redshift-space clustering with a small num¬ 
ber of parameters. The coherent bulk motion of galaxies induced by 
gravitational evolution squashes the distribution of galaxies along 
the line-of-sight direction (Kaiser 1987; Hamilton 1992; Peacock 
et al. 2001). The anisotropy in the clustering provides a good probe 
of cosmic growth rate and have been used for testing gravity (e.g, 
Guzzo et al. 2008; Yamamoto et al. 2008, 2010; Beutler et al. 2014). 
Very recently Hearin (2015) points out the impact of the assembly 
bias using color abundance matching technique and shows that the 
simple formalism of halo occupation distribution causes a signifi¬ 
cant systematics in RSD studies using pairwise velocity statistics. 

The paper is summarized as follows: in Section 2, we explain 
the details of the observational catalogs and subhalo samples used 
in our analysis. We also make a brief summary of the two color 
assignment schemes based on the subhalo age and local dark mat¬ 
ter density. In Section 3, we show the results of the luminosity and 
color dependence of observed redshift-space clustering compared 
with those of subhalos, and test SHAM and two color assignment 
methods. We discuss what causes the difference of the two color 
assignment methods. Section 4 is devoted to summary and conclu¬ 
sions. 

2 OBSERVATION AND SUBHALO SAMPLES 

2.1 Observational data 

We use the magnitude-limited samples of the Sloan Digital Sky 
Survey (SDSS) Data Release 7 galaxies (Abazajian et al. 2009) in 
three different luminosity bins: —22 < Mr < —21, —21 < Mr < 
—20, and —20 < Mr < —19 where Mr is the r-band absolute 
magnitude including —5 log h term. The apparent magnitude rrir 
ranges from 14.5 to 17.6 based on the Petrosian magnitude and K- 
corrected to rest-frame magnitude al z = 0.1. The galaxy color is 
divided into red and blue on the K-corrected g — r color and r-band 
apparent magnitude plane: red when g — r > 0.21 — O.OSMr, oth¬ 
erwise blue (see the details of the samples in Zehavi et al. 2011). 
The fractions of red galaxies become 0.63, 0.53, and 0.43 from 
luminous to faint magnitude samples. Due to the physical size of 
the spectroscopic fiber, both spectra of adjacent galaxies with the 
angular separation less than 55 arcmin cannot be measured simul¬ 
taneously, i.e., “fiber collision effect”. This corresponds to the pro¬ 
jected comoving scale Vp = 0.12 /i“^Mpc at the outer edge of the 
sample. In our analysis, we focus on the range of scales where the 
fiber collision effect is unimportant. 


2.2 Subhalo catalogs and color assignment 

We use the subhalo samples corresponding to the magnitude ranges 
of the observational samples. Below we summarize the catalogs 
briefly (see the details in Masaki et al. 2013b). The subhalo cata¬ 
logs are constructed from 77-body simulations using publicly avail¬ 
able code Gadget- 2 (Springel et al. 2001b; Springel 2005). The 
initial condition is a random-Gaussian field with the power spec¬ 
trum based on a flat A cold dark matter model of the Wilkinson Mi¬ 
crowave Anisotropy Probe (WMAPj 7-year results (Komatsu et al. 
2011): Dm = 0.272, Db = 0.0441, Da = 0.728, ITo = lOOh = 

70.2 km s“^Mpc“^,(J8 = 0.807, and Us = 0.961. We employ 
1024® dark matter (DM) particles in each cubic simulation box. We 
use two simulation boxes with the side length Lbox and the redshift 
of output snapshot Zout depending on the range of magnitude as 
listed in Table 1. 

Halos are identified using the Friends-of-Friends (FOF) algo¬ 
rithm with the linking length of 0.2 times the mean interparticle dis¬ 
tance. Satellite subhalos, dense clumps within each halo, are iden¬ 
tified using the SUBFIND algorithm with the minimum number of 
DM particles set to be 20 (Springel et al. 2001a). A central subhalo 
is defined as the rest of DM particles without satellite subhalos (and 
the ‘fuzz’ component). The central positions of both central and 
satellite subhalos are the density maximum rather than the center 
of mass, and then the velocity of a subhalo is given by the mean ve¬ 
locity of all inner particles. We neglect the inner velocity of central 
galaxies against the halo bulk velocity. 

We assume the monotonic relation between the galaxy lumi¬ 
nosity measured with Mr and the maximum circular velocity UmS 
to match their abundance as ngai(> L) = nsubhaio(> Umax) 
(Conroy et al. 2006). Umax is the maximum value of the circu¬ 
lar velocity of particles Vcirc{R) = \/ GM{< R)/R, where R 
is the distance from the center of each subhalo. The maximum cir¬ 
cular velocity for a central subhalo is computed at the observed 
epoch, however, Umax for a satellite is computed at the accretion 
epoch because subhalos inside clusters lose their initial mass due 
to tidal stripping while the stellar mass is tightly bound (Nagai 
& Kravtsov 2005). Furthermore we divide each magnitude sample 
into red and blue colors based on the following two models (Masaki 
et al. 2013b). 

(i) Subhalo age model (“Age model”) 

The galaxy color reflects the age of the stellar population: red 
galaxies consist of old stellar population when the star formation 
is quenched, whereas blue galaxies consist of younger stellar pop¬ 
ulation with active star formation. We assign the color by the sub¬ 
halo formation epoch ztoim, which is defined as the epoch when 
the maximum circular velocity of the subhalo first crosses some 
constant value / times Umax^ 

Umax(« = Zform) = / X Umax (0 < / < 1), (1) 

(ii) Local DM density model (“L.D. model”) 

The other way of color assignment is based on the local DM density 
motivated by the observation that redder galaxies locate in denser 
environments (Zehavi et al. 2005; Mandelbaum et al. 2006). The 
local DM density p{RnM) of each subhalo is estimated from the 
number of DM particles within the distance of i?DM from the sub¬ 
halo center. 

We assign color information with subhalos by assuming a 
monotonic relation with Zform, or p(i?DM) to match the red fraction 
of each observational sample. Followed by Masaki et al. (2013b), 
a constant factor / in the age model is set to be 0.9 for all magni- 
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Table 1. Basic information of subhalo catalogs for three magnitude samples: simulation box size Lpoj,. the output redshift of simulation snapshot ^out, the 
number of subhalos matching to the observed number density, satellite fraction /sat of red and blue subhalo catalogs in two color assignment models: 

subhalo age (Age model) and local DM density (L.D. model). 


Mr 

^out 

-^box 

^Buh 

All 

/sat 

Age model 
red blue 

L.D. 

red 

model 

blue 

[-22,-21] 

0.1 

300 L-iMpc 

29969 

0.22 

0.28 

0.12 

0.26 

0.16 

[-21,-20] 

0 

200 L-iMpc 

43200 

0.32 

0.40 

0.21 

0.42 

0.20 

[-20,-19] 

0 

200 L-iMpc 

81600 

0.33 

0.42 

0.26 

0.59 

0.13 


tude bins to agree with the projected clustering for both red and 
blue galaxies (Zehavi et al. 2011). The scale 7?dm in the L.D. 
model is set to be 250 fe“^kpc for the most luminous sample of 
—22 < Mr < — 21 and 200 /t“^kpc for the other samples. 

The fraction of satellite galaxies play a significant role in 
the small-scale clustering because of their large internal motion 
compared to the central galaxies. Table 1 lists the satellite frac¬ 
tion /sat = Nsa.t/{Ncen -f Nsa.t) of red and blue samples in the 
two color assignment models. In both models, red galaxy samples 
have a larger /sat than the blue sample. Luminous galaxy samples 
have smaller /sat, which reflect that the luminous galaxy is more 
likely to be central galaxies rather than satellites. The difference of 
/sat among the two color assignment models is small expect for 
the faintest sample where (smaller) satellite fraction of red (blue) 
galaxies in the L.D. model. 


3 RESULTS 

In this section, we show the result of the redshift-space clustering 
before and after color separation to see how SHAM and two color 
assignment models reproduce the observations. We use multipole 
correlation functions ^i{l = 0,2,4) and focus on the non-linear 
gravitational scales from subMpc to ~ 10 Mpc scales. 


3.1 Multipole correlation functions 

We compute the redshift-space two-point correlation functions of 
observational galaxy samples using Landy & Szalay estimator 
(Landy & Szalay 1993) 

,, , DD — 2DRRR 

=^ -, ( 2 ) 

where DD, DR, and RR is the number of data-data, data-random, 
and random-random pairs normalized by the total number of pairs 
at each bin of the three-dimensional distance of s and the cosine 
of the angle between the separation and the line-of-sight direction 
fi. Our bins of s is logarithmically equal from 0.32 to 40 /i“^Mpc 
with 11 bins. 

The anisotropy of galaxy clustering in redshift space is de¬ 
scribed with the multipole expansion using Legendre polynomials: 

Ci(s) = ^y (3) 

where jj, is the cosine of the angle to the line-of-sight direction and 
Ci{d) is the Lth Legendre polynomials: €-o{u) = 1, ^ 2 (/r) = 
(3^2 - l)/2, and Cffi) = (35Ai'‘ - 30fd + 3)/8. 


In our analysis, we focus on the monopole (I = 0 ), quadrupole 
[1 — 2), and hexadecapole (I = 4) components. Coherent bulk mo¬ 
tion of galaxies known as Kaiser effect squashes the redshift-space 
galaxy clustering along the line-of-sight direction and then the 
quadrupole component becomes negative. Kaiser effect on higher- 
order multipoles is small. On the other hand, the elongated distribu¬ 
tion by the FoG effect generates high-order anisotropic components 
such as a hexadecapole and then both quadrupole and hexadecapole 
components have positive values (Hikage & Yamamoto 2013). As 
the FoG effect increases, their amplitude becomes larger. 

Since both observational and simulated subhalo samples are 
one realization in each magnitude bin, we estimate the error of 
using the Jackknife re-sampling method as follows: 

N 

= (4) 

i 

where N is the number of subsamples, Q denotes l-lh multipole for 
the sample without data in i-th sub-volume, and fi{s)) represents 
the ensemble average of Q{s) over subsamples. In order to estimate 
the error, we divide each observational sample into 125 subsamples 
with equal sky area. Subhalo samples are also divided into 125(= 
5®) sub-cubes with equal volume. 

3.2 Luminosity dependence 

First we show the comparison of redshift-space clustering before 
color separation between the observed galaxies and the SHAM- 
based subhalo samples. Fig. 1 shows the results of = 

0, 2, 4) for three magnitude samples. The errors of both the obser¬ 
vations and the subhalo samples are estimated from the Jackknife 
resampling method. We find that the SHAM method well repro¬ 
duces the luminosity dependence of redshift-space galaxy cluster¬ 
ing without introducing additional parameters. The overall ampli¬ 
tude of increases for more luminous samples because luminous 
galaxies are hosted by massive halos and has a larger galaxy bias¬ 
ing. One can find that ^2 and ^4 are suppressed on subMpc scales in 
the most luminous sample. This is because more luminous galax¬ 
ies have a lower /sat (see Table 1) and thereby the FoG anisotropy 
due to satellites becomes smaller. The agreement of between the 
observations and the subhalos is particularly excellent for the in¬ 
termediate magnitude sample. This indicates that the internal mo¬ 
tion of satellites responsible for the FoG effect is well described 
by the motion of subhalos and then the SHAM approach can be 
a powerful tool for building a theoretical framework of redshift- 
space clustering. There is however a systematic difference that the 
brightest (faintest) galaxy sample has weaker (stronger) FoG ef¬ 
fect than the subhalo sample. This may indicate the importance 
of baryonic physics such as hydrodynamic drag by the gas that 
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Figure 1. Compaiison of the multipole correlation functions s^i (s) {I = 0,2,4 from left to right) between observations (symbols) and subhalo catalogs (lines) 
for three magnitude samples. The error-bars represent the Icr dispersion estimated with the Jackknife resampling method. Multiplication by the radial scale s 
is for a visual purpose. 



Figure 2. Color difference of the two dimensional correlation function ^(rp,7r) where Vp is the projected separation and tt is the line-of-sight separation. 
Upper (Lower) panels show the results of red (blue) subhalos by separating color by the subhalo age (left) or the local DM density (center). Right panels show 
the observational results for comparison. All figures correspond to the sample with the intermediate magnitude range of —21 < Mr < —20. 
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Age model 



I .I . . . I. I . I . . . I. I . I 
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Figure 3. Comparison of the observed multipole correlation functions ^i{s) {I = 0, 2, 4) for red and blue galaxies with those of simulated subhalo samples 
where the color is assigned by the subhalo age. Observational results for red and blue galaxies are plotted with filled circles and open triangles respectively. 
The simulation results are represented by the shade regions. Both errors of the observations and the simulations are estimated with the Jackknife resampling 
method. 


Local density model 
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Figure 4. Same as Fig. 3 but for the color assignment based on the local DM density model. 
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is not included in A^-body simulations. Hydrodynamical simula¬ 
tions show that baryonic components resists the tidal disruption of 
slowly-rotating subhalos and then decreases the averaged velocity 
of satellite subhalos (Faltenbacher et al. 2005; Wu et al. 2013). 
Our results may depend on the algorithms of finding subhalos. It 
would be interesting to apply other halo finders (e.q., Behroozi et al. 
2013), however, these works are beyond the scope of our paper and 
we leave them to future works. 


3.3 Color dependence: age model vs local DM density model 

Next we see the color dependence of redshift-space clustering and 
test two color assignment models based on the subhalo age and 
the local DM density. Fig. 2 shows the two-dimensional correla¬ 
tion functions ^{vp, tt) where the distance of galaxy pairs are sepa¬ 
rated into the tangential separation rp and the line-of-sight separa¬ 
tion TT. Upper (lower) panels show the results of red (blue) galaxies 
in the age model (left), the L.D. model (center), and observations 
(right) for the intermediate magnitude sample. Kaiser effect due to 
the coherent infalling motion squash the galaxy distribution along 
the line-of-sight direction on large scale. On small scale around 
l/i“^Mpc, the FoG effect due to the internal motion of galaxies 
inside clusters elongate the distribution of galaxies along the line- 
of-sight direction. One can clearly see the color difference of the 
FoG effect in both color assigning models: red galaxies has much 
stronger FoG effect than blue galaxies. This feature is consistent 
with the observational results for the SDSS galaxy samples in the 
right panels and also with the previous work done by Zehavi et al. 
(2011). One of the reason is that red galaxy samples have a larger 
satellite fraction than blue ones as shown in Table 1. The internal 
velocity dispersion of red galaxies is also larger because red galax¬ 
ies are hosted by massive halos. We make further discussion in the 
next subsection 3.4. 

We make a more detailed comparison using a multipole ex¬ 
pansion. Fig. 3 shows the measurements of monopole {I = 0), 
quadrupole {I = 2), and hexadecapole {I = 4) components for red 
and blue galaxies of the three magnitude samples. We find that the 
red galaxies have larger amplitudes than blue galaxies for all of the 
multipoles in our focused range of scales where the FoG effect is 
important. For comparison, we plot the results of subhalo samples 
with color assigned by the subhalo age. We find that the age model 
reproduces the color difference of redshift-space clustering down to 
subMpc scales very well. There are some systematic deviations: 
for blue subhalos in the most luminous bin are stronger, but for 
red subhalos in the faintest bin are weaker. The systematic trend is 
also seen before color separation (see Fig. 1 for comparison). This 
indicates that the discrepancy is mainly due to the incompleteness 
of SHAM not the color assignment scheme based on the subhalo 
age. We also compare the observations with the L.D. model in Fig. 
4. In contrast to the age model, the L.D. model has significant de¬ 
viations for all of the magnitude samples and multipoles at less 
than a few Mpc scale: red (blue) subhalos have too large (small) 
amplitude compared to the observations. This means that the color 
assignment by the local DM density is not enough to describe the 
redshift-space clustering and that the subhalo age is a much better 
proxy of galaxy color. 

We evaluate the agreement of the multipole correlation func¬ 
tions between the observations and the subhalo samples in the chi- 
squared basis. The total in each magnitude bin is computed by 
summing up the x^ for red and blue galaxies over different multi¬ 
poles Z = 0, 2, and 4: 


Table 2. Reduced chi-squared values of the multipole power spectra (Z = 
0, 2, 4) for three magnitude samples before color separation (“No color 
sep.”) and two color models (Age and L.D. model). Each chi-squared value 
is summed over three multipoles I = 0,2, and 4 and two colors after 
color separation. We adopt the range of scales is from 0.32 Zi~^Mpc to 
40 /i“^Mpc and the number of bins are 11 in the range. The total d.o.f. is 
33 (or 66) before (or after) color separation. 


Mr 


X^/d.o.f. 


No color sep. 

Age model 

L.D. model 

[-22,-21] 

3.49 

3.52 

8.45 

[-21,-20] 

0.89 

1.87 

6.11 

[-20,-19] 

4.01 

2.77 

6.18 

8^- ^^^^^^^- 


o 

■a 


<1 


Local density • 
Age -A- 


-A-. 


Ax 


s[h’^Mpc] 


Figure 5. Comparison of the reduced chi-squared values between the age 
model and the L.D. model on different bins of scale s. Lower panel shows 
the difference of the chi-square Ax^ between the two models. 


0,2,4 <s j <[Smax 

2 

X = 


color I 






(5) 


1 obs, 2 1 model, 2 .1 • r /- r 1 

where cr^ ’ , and cr^ are the variance or tor observa- 

tions and subhalo samples estimated from the Jackknife method. 
The range of the fitting scale is from Smin = 0.32/i“^Mpc to 
Smax = 40Zi“^Mpc with 11 bins . We neglect the covariance be¬ 
tween different scales, multipoles and colors for simplicity. 

Table 2 lists the reduced chi-squared values of three magnitude 
samples before and after color separation. The degrees-of-freedom 
(d.o.f.) of each magnitude bin is 33 (11 bins of scale x 3 mul¬ 
tipoles) before color separation and 66 after color separation. We 
find that the reduced x^ of the age model is 2-3. On the contrary, 
the L.D. model is significantly larger chi-squared values such as 
6-8. Fig. 5 shows the comparison of the reduced x^ at each bin of 
scale on l/i“^Mpc between the two color models. The fitting 
of the L.D. model is already worse around l/i“^Mpc and then the 
deviation increases as the scale goes down. 


3.4 Discussion 

In the previous subsection, we find that the age model has a bet¬ 
ter agreement with the observed redshift-space clustering than the 
L.D. model. We here investigate why the two models have different 
features in the redshift-space clustering. The amplitude of from 
subMpc to Mpc scale is mainly affected by the FoG effect, which 
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Figure 6. One-dimensional internal velocity dispersion (tid of red (left) and blue (right) satellite subhalos relative to the central subhalos as a function of 
the host halo mass for three magnitude samples. The color separation is based on the age (upper) and local DM density (lower) models respectively. For 
comparison, Virial velocity dispersion is plotted with lines. The eiTor-bar represents Icr dispersion in each bin of halo mass. The lower panels show the ratio 

O’lo/o'vir- 


depends on the fraction of satellite galaxies, the internal velocity 
dispersion, and also the radial profile of satellites in their host ha¬ 
los. We here investigate where the differences between the models 
come from one by one. 

As the satellite fraction increases, the contribution of central- 
satellite pairs and satellite-satellite pairs in the same halos in¬ 
creases, which generate the large amplitude of quadrupole and hex- 
adecapole components on ^ IMpc. As shown in Table 1, however, 
the satellite fraction between the two color models is roughly same 
except for the faintest magnitude bin and thereby the number of 
satellites is not a main reason for the model difference of redshift- 
space clustering. 

The FoG effect comes from the internal velocity dispersion, 
which is mainly dependent on the host halo mass. Table 3 lists the 
mean halo mass weighted with the number of satellite subhalos: 


/ Ti^sat,wei\ _ 

V^'^halo / ~ 


-^^halo,i-^sat,i 


( 6 ) 


where Mhaio,i and Asat.i is the mass and the number of satellite 
subhalos in i-th halo. One can see that more luminous and redder 
satellites are hosted by more massive halos. We find that the differ¬ 
ence of the host halo mass between colors is rather small in the L.D. 
model than the age model. It is thereby difficult to explain the large 
color difference of redshift-space clustering in the L.D. model. 

We also see if the velocity dispersion of red and blue sub¬ 
halos follow the expectations of Virial theorem. Fig. 6 shows the 


Table 3. Comparison of the mean halo mass weighted with the number of 
satellite subhalos (eq.[6]) for the red and blue subhalo samples in the two 
color assignment models. 


Mr 


[lOi-^/i-iMo] 

Age model 

L.D. model 

red 

blue 

red blue 

[-22,-21] 

1.46 

1.09 

1.40 1.36 

[-21, -20] 

1.30 

0.80 

1.13 1.16 

[-20, -19] 

1.06 

0.64 

0.90 0.76 


one-dimensional (ID) internal velocity dispersion of satellite sub¬ 
halos as a function of the host halo mass in the age and L.D. 
model. Lines represent the Virial velocity dispersion estimated by 
cTvir = (GiMvir/Si^vir)^^^ where Mvir and i?vir represent Virial 
mass and radius. The ID velocity dispersion of satellite subhalos 
to the central subhalos are estimated as The 

error-bars represent the 1 -ct dispersion in each bin of halo mass. 
In both models, we find that the satellite velocity dispersions agree 
with the expectations from Virial theorem within the error bars ir¬ 
respective of color, though dispersion slightly decreases as the host 
halo is more massive. The velocity dispersion is unlikely to explain 
the difference of redshift-space clustering between the two models. 
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Figure 7. The fraction of red galaxies as a function of the distance from the host halo center normalized by the Virial radius iJvir for three magnitude samples. 
Red and blue lines represent the results for the L.D. model and the age model respectively. The error-bar represents the Poisson error in each bin of scale. 


Finally we see the radial profile of satellites which affect the 
small-scale clustering. Fig. 7 shows that the fraction of red galax¬ 
ies /red = Nred/iNved + Nhiue) as a function of the distance R 
from the halo center normalized by the Virial radius J?vir in the 
three magnitude samples. For all of the samples, there is a clear 
difference between the two color models. The red fraction in the 
age model slightly decreases at larger R, whereas the red fraction 
in the L.D. model is almost unity inside the host halo and dras¬ 
tically declines around the Virial radius. This can explain the big 
difference between colors in the redshift-space clustering: a large 
number of red galaxy pairs inside halos causes too strong FoG ef¬ 
fect, however, the lack of blue galaxies have too weak FoG. The ob¬ 
served red fraction actually shows a similar trend to the age model 
(Hansen et al. 2009). This means that the color assignment in the 
L.D. model is too simple to describe the color difference of galaxy 
clustering inside halos. 

We investigate how the result changes by including scatter in 
the color assignment. Instead of assuming the monotonic relation 
between the color and pdm, we assign red color to the subhalo with 
the local DM density of pdm in the following probability: 


f’red(pDM) = 


1 + erf 


^^°glo(PDM/pD*M ) 

Criog p 


(7) 


where erf is the error function, Pi^M is a threshold density divid¬ 
ing red and blue colors without scatter to match the observed red 
fraction, and criogp denotes the scatter in logarithmic scale of pdm- 
Fig. 8 shows the comparison of observed quadrupole components 
of the intermediate sample with the corresponding measurements 
in the L.D. model by varying uiog p as 0.2, 0.3, and 0.4. We change 
the parameter of i?DM to be 600 kpc to agree with the observed 
projected correlation functions (Masaki et al. 2013b). Introducing 
an appropriate scatter improves the agreement with the observa¬ 
tions. We find that the scatter value of criogp around 0.3 has a best 
agreement in blue samples, while criogp between 0.3 and 0.4 is best 
for red samples. The chi-squared value for the criogp = 0.3(0.4) 
significantly decreases to be 3.62(4.55) compared to the no-scatter 
case of 6.11(6.18) in Table 2. The agreement in the age model is 
still better than the L.D. model even when the scatter is included. 
This is because the introduction of scatter in the form of eq. [7] 
cannot reproduce the clustering of red and blue galaxy samples si¬ 
multaneously. This implies that the age model is much better color 
assignment than the L.D. model and then the subhalo age is a key 
ingredient to determine the color of galaxies. 


4 SUMMARY AND CONCLUSIONS 

We extend the SHAM approach to redshift space and test if 
the SHAM explain the luminosity and color dependence of the 
redshift-space clustering. We find that the simple subhalo abun¬ 
dance matching using monotonic relation between the galaxy lu¬ 
minosity to the maximum circular velocity qualitatively well repro¬ 
duces the luminosity dependence of SDSS galaxy clustering from 
0.3 /i“^Mpc to 40 /i“^Mpc. This indicates that the satellite motion 
inside clusters is mainly determined with that of the host subhalos. 
Our results indicate that the SHAM method can be applied for RSD 
studies and provides a promising way to construct mock samples of 
redshift-space galaxy distribution. There is however a systematic 
difference on subMpc scale and then the effect of baryonic physics 
may be necessary to be included to achieve more precise theoretical 
modeling for future galaxy surveys. 

In addition to the luminosity dependence, we also apply two 
methods of color abundance matching where the color is assigned 
by the age of subhalos and the local DM density. We find that the 
color assignment by the subhalo age much better agrees with the 
observations than by the local DM density. The main reason why 
the local density model fails to reproduce the observed clustering 
is that the fraction of red subhalos in the model is too large inside 
the host halos. The agreement improves by introducing a scatter in 
the relation between color and local DM density, but still the color 
assignment based on the subhalo age has better chi-squared values. 
This suggests that the subhalo age is a main driver of determining 
the color of galaxies and a key ingredient to make an accurate mock 
galaxy samples with different colors. 

Our analysis using redshift-space clustering prefers the age 
model to the local density model. This is apparently inconsistent 
with the galaxy-galaxy lensing analysis by Masaki et al. (2013b), 
which support the local DM density model. This is because the 
redshift-space clustering and lensing are sensitive to different as¬ 
pects of the galaxy clustering: the redshift-space clustering, in par¬ 
ticular the FoG effect, is sensitive to the satellite properties such as 
the satellite profile and dynamics. On the other hand, galaxy-galaxy 
lensing probes the averaged halo mass hosting all of the galaxies in 
a given sample. This means that both of the color assignment mod¬ 
els still need to be improved. 

In this analysis, we assume that the central galaxy sit on the 
halo center and neglect the internal velocity relative to the host 
halo bulk velocity. Recent analysis using the Baryonic Oscillation 
Spectroscopic Survey (BOSS) CMASS data shows that the cen¬ 
tral galaxies may have ~ 30% of the Virial velocity (Guo et al. 
2015). This increases the velocity dispersion between central galax¬ 
ies and satellite galaxies at ~ 9% (~ 0.3^). Baryon components 
also affect the dynamics of satellite galaxies. The mean velocity of 
satellites decrease by including baryonic components which pre- 
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Figure 8. Compaiison of observed quadrupole component ^2 for red and blue galaxies based on the L.D. model. The scatter in the color assignment is included 
as the equation 7 and the logarithmic scatter value triog p is 0.2, 0.3, and 0.4 from left to right panels. 


vents the tidal disruption of slow subhalos by 10 percent level (Fal- 
tenbacher et al. 2005; Wu et al. 2013). It may be also interesting to 
see how our result changes using different algorithms of identifying 
subhalos. We leave further detailed analysis as a future work. 
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